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We give an update on our ongoing investigations of potentials between pairs of static-light 
mesons, ^(r)^(O) and ,?^(r),?^(0), in A^f = 2 Lattice QCD, in different spin and isospin channels. 
The question of attraction and repulsion is particularly interesting with respect to the X(3872) 
charmonium state and charged candidates such as the Z+(4430). We employ the nonpertur- 
batively improved Sheikholeslami-Wohlert fermion and the Wilson gauge actions at two lattice 
spacings a « 0.084 fm and a « 0.077 fm with a pseudoscalar mass of mps ~ 770 MeV and 
mps « 400 MeV respectively. We use stochastic all-to-all propagator techniques, improved by a 
hopping parameter expansion. The analysis is based on the variational method, utilizing various 
source and sink interpolators. 
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1. Introduction 

Potentials between static-light mesons = Qq) are of interest since they give insights in the 
nature of strong interactions from first principles for multiquark systems. For large heavy quark 
masses, e.g., the spectra of heavy-light mesons are determined by excitations of the light quark and 
gluonic degrees of freedom. In particular, the vector-pseudoscalar splitting vanishes and the static- 
light meson 3S can be interpreted as either a B, a B*, a D or a D* heavy-light meson. Calculating 
potentials between two ^ mesons then will also enable investigations of possible bound tetraquark 
states or for particles that are close to the meson-antimeson thershold, such as the X(3872) or 
the Z+(4430). Therefore, and .'^^E potentials have been studied by many groups. First 
calculations were performed by Michael and Pennanen [1,2]. A more detailed quenched study 
can be found in ref. [3]. For the computation of static-light meson-antimeson systems we refer 
to [4]. Recent dynamical simulations with twisted mass fermions were carried out by Wagner in 
refs. [5, 6] and with Sheikholeslami-Wohlert fermions in our Lattice 2010 proceedings [7]. 



2. Computation 

In these proceedings we present the latest results from our investigations of potentials between 
two static-light mesons. Of special interest is the question of attraction and repulsion and their 
dependence on the sepai^ation between the static quarks. Therefore, we numerically determine 
ground and excited states of ^ mesons as well as intermeson potentials between pairs of static- 
light mesons, =^(r)=^(0) and <^(r)=^(0). The static quark-quark (or quark-antiquark) separation 
is given by r = |r| = Ra,R G INq. a denotes the lattice spacing, Q a static colour source and the po- 
sitions of the mass-degenerate light quarks q G {u,d} are not fixed. Thus, ^ mesons carry isospin 
7=1/2 and for and states 7 G {0, 1} with 7^ G {-1,0, 1}. For =^ mesons the isosinglet 
corresponds to the representation q^qj = ud + du and the isotriplet to qiq2 G {uu,dd,ud — du} 
respectively. The corresponding representations for states are given by q\q2 = uu + dd for 
7 = and q\q2 G {du,iid, uu — dd} for 7 = 1. Graphically, the quark line diagrams with respect to 
isospin that we evaluate can be depicted as, 

Isospin 7 = 0: CMt) =^^^ Cfijwi^) ' 

Isospin7 = l: C^^(0=|j^-^, ^^^(0=^(5^- (2-2) 

Straight lines represent static quark propagators and wiggly lines light quark propagators. 

Masses are calculated from the asymptotic behavior of Euclidean-time correlation functions. 
We employ N{ = 2 Sheikholeslami-Wohlert configurations generated by the QCDSF Collabora- 
tion [8]. The parameter values are listed in table 1, where the scale is set using ro(jS, k) = 0.5 fm. 



volume x 




'Qal — 'Qea 


csw 


a/fm 


La /fm 


mps /MeV] 


^conf 


16^ x32 
24^ X 48 


5.29 
5.29 


0.13550 
0.13620 


1.9192 
1.9192 


0.084 
0.077 


1.34 
1.85 


770(9) 
400(4) 


180 
200 



Table 1: Lattice parameters. 
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The pseudoscalar mass corresponds to its infinite volume value. We use the Chroma software sys- 
tem [9]. For the techniques and improvement methods we use, we refer to [7] and the references 
therein. To analyze our data and to extract also excited states we apply the variational method [10], 
solving a generalized eigenvalue problem for a 3 x 3 cross correlation matrix generated by different 
amounts of Wuppertal smearing [1 1] applied to the source and sink operators. EiTors are calculated 
using the jackknife method. 

3. Representations and classification of states 

To create =^ meson states as well as .'^SS and systems of different y^^*-' we use interpo- 
lators ^ = Qi?q, where the operators 0" contain combinations of Dirac /-matrices and covariant 
lattice derivatives. This has been discussed in [7] and we give a summary for the representation 
and classification of our states. 

In the continuum limit, the static-light states can be classified according to fermionic representa- 
tions of the rotation group 0(3). At vanishing distance r = the and states can be 
characterized by integer and J^'~ quantum numbers, respectively. However at r = |r| > the 
0(3) (or 0(3) (X" '^) symmetry is broken down to its cylindrical Dooh subgroup. The irreducible 
representations of this are conventionally labeled by the spin along the axis A, where £,11, A refer 
to A = 0, 1 , 2, respectively, with a subscript T] = g for gerade (even) PC = + or Tj = m for ungerade 
(odd) PC = — transformation properties with respect to the midpoint. All A > 1 representations are 
two-dimensional. The one-dimensional £ representations carry an additional Ov = ± superscript 
for their reflection symmetry with respect to a plane that includes the two endpoints. 

The operators that we used to create the static-light mesons are displayed in table 2. The in- 
termeson potentials were obtained by combining two static-light mesons of different (or the same) 
quantum numbers. This can be projected into an irreducible Dooh representation, either by coupling 
the light quarks together in spinor space [5] or by projecting the static-light meson spins into the 
direction f of the static source distance, by applying ^(1 it //sy- f), and taking appropriate sym- 
metric (Aj = 1) or antisymmetric {K- = 0) spin combinations. These two approaches can be related 
to each other via a Fierz transformation. For our coarse lattice and the operator combinations that 
couple to total angular momentum / = we have performed this projection. For the other com- 
binations and our fine lattice different representations will mix. The analyzed operators and the 
corresponding representations are listed in table 3. We note that the operator combinations 75 x 75 
and 1x1 carry the same quantum numbers as well as the combinations 75 x 1 and 75 x jiVi. 





wave [12] 


Oh' rep. 


continuum 


(heavy-light) 


75 


5 


Gt 


1 + 
2 


0,1 


1 


P 


Gx 


1 - 
2 


0+,l+ 


7Vr 


P 




1 - 
2 


0+,l+ 


(7iVi -72V2)+cycl. 


P+ 


H 


3- 
2 


l+,2+ 



Table 2: Operators and representations for static-light mesons. In the last column we display the for a 
heavy-light meson, obtained by substituting the (spinless) static source by a heavy fermion. 
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Table 3: Operators and continuum representations for the meson-meson and meson-antimeson (MM) 

potentials in the isosinglet and the isovector channel. 

4. Results 

The eigenvalues X^'^^tjo) of the generalized eigenvalue problem [10], are fitted to one- and 
two-exponential ansatze, to obtain the ^h mass. The appropriate values of to and the fit ranges 
in t are determined from monitoring the effective masses as described in [7]. Let us first discuss 

meson systems. We define intermeson potentials as the differences between the meson-meson 
energy levels and the r — )■ 00 two static-light meson limiting cases: 

%ii^2 ('■) = ■£!«|.«2 ('') - ("^.^l + '".!^2 ) ^ • (4-1) 

In the left panel of figure 1, we display the ground state (E+) and the first excited state (Sg ') of the 
75 ^ 75 operator as well as the £,7 ground state and the £,7' first excited state of the 75 (g) 1 operator, 
both in the / = 1 , Aj. = channel. In the last case the lowest lying r — )• 00 combination of states 
would be a radially excited ^ state (75) plus a ^ (1) ground state. The next level would be 
the sum of and ^ It is not clear to which one of these states our creation operator has best 
overlap. In the figure we display both possibilities. The latter assignment would mean that in the 
excited state channel (like for the ground state) there is repulsion at intermediate distances. 

In the right panel of figure 1 we show the 75 (g) 75 operator in the different spin A, = 0, 1 and 
isospin 7 = 0, 1 channels for the ground state. For short distances we observe attraction in all spin 
and isospin channels. In fact at very short distances we find attraction in all analyzed channels, see 
table 3, for ground and excited states. This may not be too surprising as this is expected from gluon 
exchange in the 3* channel between the two static sources. When comparing the ground state of 
the 75 75 operator for different spin and isospin channels we figure out that the Aj = channel is 
more attractive than the A. = 1 channel for isospin 7 = 0. For isospin 7=1 this pattern is reversed. 
In both cases the difference is of the order of 150 MeV at a distance of 0.11 fm. For the other 
spin-projected operator combination (1^1, 75 (gi 1 and 75 JiVi) we also find attractive forces of 
similar sizes for the A, = 0, 1 and isospin 7 = ground states while for isospin 7=1 the A, = 1 
channel is more attractive than the Aj. = channel. 

On our coarse lattice we observe repulsive potentials at distances between 0.2 fm and 0.45 fm 
for the ground state of the 75 (8) 1 operator in all spin and isospin channels. In addition we find 
repulsion in the 7 = 1 , A^ = 1 channel for the E+ ground state of the 75 75 operator at distances 
between 0.3 fm and 0.45 fm. In the case of the fine lattice we did not perform the A^ projection so 
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Figure 1: Intermeson potentials V{r) from our coarse lattice. The combinations 75 (g) 75 (E+,E+') and 75 (g) 1 
(Z^ , for / = 1 , Aj^ = are shown on the left hand side, the operator combination 75 (S) 75 in different spin 
Aj^ = 0, 1 and isospin / = 0, 1 channels on the right hand side. 
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Figure 2: Comparison of intermeson potentials V{r) in different isospin channels between the results from 
our coarse and our fine lattice and heavy and light pion masses respectively. The combination 75 (g) 75 (E+) 
is shown on the left hand side, the combination 75 g) 1 (I.^) on the right hand side. 



that here we cannot distinguish between £ and IT states. In agreement with the coarse lattice results 
we obtain repulsion of 0(50 MeV) in the ground states of the operator combinations 75 (g) 1, 75 (gi 
V,7 (Ej^) and 75 (g) 75 (£+) for isospin 7 = 1 at intermediate distances r > 0.35 fm. 

In figure 2 we compare coarse and fine lattice results in the 75 (g 75 (£+ ground state) and 75 (g 1 
(£,7 ground state) channels. In the first channel we observe reasonable scaling while in the latter 
channel the fine lattice potentials appear to be more attractive. This may be related to the lighter 
pion mass resulting in a different Yukawa interaction. 

In figure 3 we display the 7 = 0, Aj = ground state for the meson-antimeson case in the 
75 (g 75 (£+) channel. On the left hand side we see the effective ground state energy levels E^ff for 
different ?o of the 75 (g 75 operator combination at a distance of 0. 11 fm. One finds a short plateau of 
poor quality in a range f/c? G {4, . . . , 8} for fo > 2a. Then the effective energy level decreases again 
and forms another plateau from t/a = S onwards. This can be explained by the observation that 
this state has the same quantum numbers as the QQ static potential (and multipaiticle states of the 
static potential plus a P wave rj meson, the static potential plus 2 pions etc.). Our interpolator basis 
however, has very little overlap with these states. Therefore, we cannot easily disentangle the QQ 
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Figure 3: Effective masses for the 75 ® 75 SS,'^ isoscalar Aj^ = ground state for different tQ (left hand side). 
Masses for the 75 ® 75 ground state. The horizontal lines denote the sum of two static-light meson 
masses while the black symbols are the static QQ potential plus various light meson masses (right hand 
side). 



Static potential and this background of multiparticle excitations from the lowest lying SS^^ state that 
we are interested in. The ansatz to rewrite the correlator as Cij{t) = 0^'^{t) + dij ■ exp(— • t) 
and extract Cfj^(t) failed as well. Thus, we tried to fit the ground state mass without using 
the variational approach, but from our single correlation function that has the largest overlap with 
the c^c^ ground state. Masses could be extracted for separations r < 0.26 fm. For larger distances 
the coupling to the state could not be resolved. The result is displayed in the right panel of 
figure 3 (blue circles). The two horizontal lines correspond to twice the ground state mass of the 
^ static-light meson, the expected r — )• 00 limit. At first sight there appear to be very substantial 
short distance attractive forces in this channel. Also the static potential Vqq is lying much lower 
and can be disentangled. However, states consisting of a static potential and a scalar 7 = particle 
will have the same quantum numbers. For our lattice parameters the P-wave pseudoscalar 77 meson 
(that at our quark mass will have a similar mass to the pion) is the lowest such state, with masses of 
a /o meson as well as two pseudoscalars lying higher. We include these sums in the figure where 
we approximate by m„ and m by mag . The ground state lies between these states and the 
static quark potential. So it is hard to decide whether we see a substantial attraction between the 
static-light meson-antimeson pair in this channel or bound states between the static quark potential 
and additional light mesons. 



5. Conclusions 



We investigated interactions between pairs of static-light mesons and found attraction for short 
distances in all spin and isospin channels. For distances of the order of 0.4 fm some operator 
combinations yield repulsion, in particular the combination 75 (g) 1. The interaction ranges are 
larger on the fine lattice with a smaller pion mass than on the coarse lattice. Meson-antimeson 
potentials are also very interesting with respect to charmonium threshold states [13] (DD molecules 
or tetraquarks) but difficult to disentangle from mesons that are bound to the static potential (hadro- 
quarkonium [14]). Analyzing 7 = 0, A, = states is a very challenging task since they couple 
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directly to the QQ static potential, e.g. in the case of the 75 (g) 75 combination, or to the vacuum 
state, e.g. for the 75 (8) 1 channel. 
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